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Abstract 

We investigate the gluon Green's function in the high energy hmit of QCD using a 
recently proposed iterative solution of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
equation at next-to-leading logarithmic (NLL) accuracy. To establish the applica- 
bility of this method in the NLL approximation we solve the BFKL equation as 
originally written by Fadin and Lipatov, and compare the results with previous 
studies in the leading logarithmic (LL) approximation. 



1 Introduction 



The Balitsky-Fadin-Kuraev-Lipatov [1] formalism resums a class of loga- 
ritlims dominant in the Regge limit of scattering amplitudes where the centre 
of mass energy ^/s is large and the momentum transfer is fixed. Within 
this approach the high energy cross-section for the process A + B ^ A' + B' 
can be written as 



- / Sfe / S *^(^) ^^-f^') / k.. A . in A) , (1, 
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with ^A,B being the process-dependent impact factors and / (ka, k^, A) the 
process-independent ghion Green's function. This ghion Green's function de- 
scribes the interaction between two Reggeised gluons exchanged in the t- 
channel with transverse momenta ka (,, and it carries the energy dependence 
of the cross-section. We choose to work with the symmetric Regge scale 
So = kakb, where ki = |kj|. 

The resummation of terms of the form (a^A)" defines the LL accuracy while 
the inclusion of contributions proportional to (ttgA)" leads to the NLL ap- 
proximation. The ladder structure of the scattering amplitudes derived in the 
BFKL formalism is described by an integral equation for the gluon Green's 
function. The eigenfunctions of this integral equation are known at LL ac- 
curacy and it is therefore possible to fully reconstruct the solution. One of 
the motivations to include the NLL contributions is to introduce running cou- 
pling effects. The logarithmic dependence introduced by the running coupling 
terms in the NLL approximation significantly complicates the study of the 
equation [2]. During the past few years different strategies have been sug- 
gested to study the NLL BFKL Green's function. In a fixed coupling analysis 
in Ref. [3] it was first highlighted that the NLL corrections are large and neg- 
ative compared to the LL. Different approaches to improve the convergence 
of the series expansion at NLL accuracy have been proposed [4]. When run- 
ning coupling effects are taken into account the situation improves as it has 
been shown in Ref. [5], in particular when those terms proportional to (3q are 
resummed into a^. 

Recently, we proposed [6] to use an iterative approach to solve the equation 
in the NLL approximation. A similar method was first suggested at the LL 
accuracy in Ref. [7] , where it was shown to reproduce the analytic solution, and 
opened the possibility to perform detailed LL phenomenological studies [8]. 

The NLL formalism of Ref. [6] has the advantage of dealing with the BFKL 
kernel as calculated in Ref. [9,10] with no approximations. In particular, the 
solution includes all running coupling effects and, as we do not use the angular 
averaged kernel, it allows for a complete study of angular dependences. The 
main purpose of this paper is to establish the applicability of the NLL solution 
put forward in Ref. [6]. With such intention we will solve the NLL BFKL 
equation as originally written in [9] , with a particular choice of renormalisation 
scale. Running coupling effects therefore correspond to an expansion of the 
one-loop running of as{ii) in the MS scheme. A study of different schemes for 
the running of the coupling and choices for the renormalisation scale will be 
presented in a future publication. 

This paper is organised as follows: In Section 2 we present the BFKL equa- 
tion and sketch the derivation of the iterative solution following Ref. [6]. In 
Section 3 we present a numerical analysis of the NLL kernel, briefly indi- 
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eating the mathematical expressions used in our mimerieal implementation. 
In Section 4 we use this numerical implementation to study the NLL gluon 
Green's function. We analyse the convergence of the solution and its depen- 
dence on the transverse momenta of the gluons exchanged in the t-channel. 
We present results on the evolution of the gluon Green's function and a study 
of the renormalisation scale dependence. We also show results on angular de- 
pendences and a toy cross section obtained using simplified LL impact factors, 
to finally present our conclusions in Section 5. 



2 The Solution of the NLL BFKL Equation 

In this Section we sketch the method of solution for the BFKL equation at 
NLL accuracy proposed in Ref . [6] . To write the equation in a convenient way 
we first perform a Mellin transform of the gluon Green's function, i.e. 

a+ioo 

/(k„,k,,A) = — J duje^'^UiKM)- (2) 

a— ioo 

With such a transformation the BFKL equation in dimensional regularisation 
(D = 4 + 2e) reads [9] 

{K, = S^'+'^^ {K -^b)+l d'+'^k' /C (k„, k') U (k', kb) , (3) 

where the kernel is expressed in terms of the gluon Regge trajectory and the 
real emission component in the following way 

/C (k„, k) = 2 cu(^) (k^) 5(2+2^) (k„ - k) + ICi^^ (k„, k) + ICr (k„, k) . (4) 

Both e-dependent parts of the kernel contain 1/e and l/e^ poles. The cancel- 
lation of the poles in the trajectory against those in the real emission kernel 
was shown in Ref. [6] by introducing a phase space slicing parameter A. The 
BFKL equation can then be expressed as 

- ouo {kl, X')) U (k„, k,) = (5(2) (k„ - k,) (5) 
+ / dF^^ (k') iy? - A^) + ICr (k„, k, + k)) (k„ + k, k,) , 

where /C^ is the finite part of the emission kernel. 
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(7) 



with q;s(//) = as{iJ,)Nc/7r, plays a crucial role in the cancellation of infrared 
singularities in the final result. It contains the information about the running 
coupling effects encoded in the terms proportional to /?o = — |^/- This ^ 
function can be modified to resum these effects and consider them in different 
schemes. We will study this in a separate publication. In the present work we 
take ^ as in Eq. (7), which corresponds to the one-loop expansion of q;s(//) in 
the MS scheme with renormalisation scale /j,. In this approach the expression 
for the gluon Regge trajectory reads 



(q^A 



2^--e(|q|A)ln^ + a,2(/z)-C(3). 



(8) 



As explained in Ref. [6] one can solve the BFKL equation by iterating Eq. (5) 
in such a way that an extra cj-pole is generated per rung in the BFKL ladder. 
It is then possible to perform the inverse Mellin transform to find the solution 
directly in energy space, i.e. 



/(k„, k,, A) = exp {uo {kl, X') a) I S^'\K - K) 



(9) 
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where we have used the notation yo = A. 

In the numerical implementation discussed in this work we have chosen the 
renormalisation scale to be // = kb, one of the perturbative scales in the 
interaction. The study of alternatives to this choice will be presented elsewhere. 

It is important to realise that Eq. (9) gives the correct solution in the A — > 
limit. In practice we are able to numerically check the region of stability of 
the expression in Eq. (9), i.e. the region at small A where the result is flat 
in A. Every extra term in the series expansion corresponds to an additional 
iteration of the kernel in the integral equation. For a given value of the variable 
A and the slicing parameter A only a finite number of terms in the expansion 
is needed to obtain the solution to a given accuracy. 

In the next Section we show the mathematical expressions used for the kernel 
in our implementation, and briefly describe the structure of the trajectory and 
the real emission kernel. 



3 Analysis of the BFKL Trajectory and Kernel 



As we have already pointed out, there are two key elements in the BFKL equa- 
tion: the gluon Regge trajectory, ujq (q^, A^), which contains the function C,, and 
the finite real emission component of the kernel, ICr (q, q'). The trajectory is 
not infrared-finite and therefore carries a A dependence. It is interesting to 
compare the behaviour of the trajectory at LL to that at NLL. This is shown in 
Fig. 1, where first we display the dependence on A with q = 20 GeV and then, 
for a fixed value of A = 1 GcV, wc plot the behaviour in q. The values used 
in these plots are n/ = 4, Ac = 3, and, for the top four graphs, — 30 GeV. 
For a value of A^^^ of 0.1416 GeV this imphes that as{n) = 0.1408. For this 
value of the renormalisation scale the NLL trajectory is always more negative 
than the LL one with the trajectories separating from each other as the slic- 
ing parameter decreases for a fixed q, or, when we fix A, for large values of 
q. For fixed values of q and A, and for large renormalisation scales the effect 
of changing does not affect the difference between the trajectory at LL and 
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NLL, as can be seen in the bottom two plots in Fig 1. This is not the case for 
lower values of fi, where the trajectories can even overlap at the minimum of 

, ,NLL 
UJq 

To implement the solution for the gluon Green's function numerically, the 
kernel in Eq. (6) is rewritten as 
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where 9 is the angle between the two-dimensional vectors q and q', and 



T{q,q\e) = Im|4Li2||^e- 



— In —n In 



q^ q' \sm.6\ — i {q — q' cos( 



q' 5' |sin^| + i (g — g' cos6') 



(11) 



with 



LM.) = -/aMLz^). (12) 



This kernel has an integrable singularity in the q — q' limit. As an example 
of its structure in the vicinity of this singularity we plot ICr {q, q' = 20 GeV, 9) 
in Fig. 2 with the same values for the parameters as those in Fig. 1 where the 
renormalisation scale was fixed to /i — 30 GeV. 

In the next Section we use this form of the kernel in order to study the solution 
to the NLL BFKL equation. 



4 Study of the Gluon's Green Function 



We have implemented Eq. (9) in a Monte Carlo integration routine using the 
form of the kernel presented in Section 3. In the following analysis, we have 
chosen to run the coupling at one loop, matching the values of Ref. [11] and, 
as previously explained, we use = /c^ as the renormalisation scale and set 
Uf — A. The function ^ will always correspond to Eq. (7). 
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Fig. 2. Structure of the kernel Kr {q, q', 6) (in GeV ^) for q' = 20 GeV as a function 
of q and the angle between q and q'. 



Convergence of the Solution 



Before presenting the results for the gluon Green's function, we first investigate 
the properties of convergence for Eq. (9). There are two points of interest: 
Firstly, how the A — >^ limit is approached, and secondly, to determine how 
many terms in the infinite sum of Eq. (9) are needed to obtain the solution 
within a given numerical accuracy. These two points are linked because the 
smaller the value of A, the more terms are necessary in the expansion to 
reach good accuracy, but also the better the approximation /^^ (k^ + k, k;,) ~ 
/t^(ka,k6) for |k| < A (as shown in Ref. [6] this approximation is used in 
order to write the BFKL equation as in Eq. (5)). A good choice for A is 
therefore characterised by being small enough for the approximation to be 
valid, ensuring in this way an accurate answer, while being large enough to 
warrant the contribution from only a finite number of terms in the expansion, 
so that the numerical evaluation is fast. 

In Fig. 3 we plot the contribution from successive terms in the infinite sum of 
Eq. (9) to the angular integrated NLL gluon Green's function 
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Fig. 3. Distribution on the number of iterations in building up the NLL gluon 
Green's function for different values of the parameter A, at a fixed value of 
A = 1 GeV. The gluon Green's function is evaluated for ka = 25 GeV, kh = 30 GeV 
and the renormalisation scale is chosen to he fi = ki,. 



with 6 being the angle between and k^, at ka = 25 GeV, kb = 30 GeV, 
A = 1 GeV, and for different values of the parameter A. These values have 
been chosen with no other intention but to illustrate the capability of the NLL 
formalism proposed in Ref. [6] and its numerical implementation. We see that 
for a given choice of parameters, only a finite number of terms contribute to the 
infinite sum. All the results presented in this paper have been calculated with 
some upper limit on the number of terms included in the infinite sum of Eq. (9). 
It has been verified that this limit is put sufficiently high as to reproduce the 
solution with the required accuracy. This plot contains information about how 
the emission builds up, revealing in a quantitative way the fact that when the 
energy available for the scattering process is larger, the distribution peaks at 
larger values of the number of iterations. Although this trend is independent 
of A, the specific position of the peak is not. 

f{ka, kb, A) in Eq. (13) could also have been obtained by first angular averaging 
the BFKL kernel and then solving the equation. Although in our solution we 
have all the angular information, in this case we average in angles in order to 
compare with the analytic expression for the LL BFKL gluon Green's function 
with zero conformal spin, i.e. 




(13) 
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(14) 



with the LL eigenvalue being 




(15) 



The LL hmit of Eq. (9) is given by 
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In this study we have checked that the LL analytic results from Eq. (14) 
coincide with those from the LL version of our numerical implementation. 
This will be illustrated in the plots below. 

Having shown how it can be checked that a sufficient number of iterations in 
Eq. (9) has been performed for a given choice of the parameters, we will now 
proceed to study the A — limit of the equation. In Fig. 4 we have plotted 
the A-dependence of the angular integrated gluon Green's function on A for a 
choice of ka — 25 GeV, kt, = 30 GeV and A = 3. The result is very flat in A for 
small values of this parameter, demonstrating, remarkably, the cancellation 
between the infrared divergences present in the gluon NLL Regge trajectory 
and those stemming from the integration of the real NLL emission kernel 
over phase space. For larger values of A we observe a growing A-dcpcndence 
originating from our initial approximation /^^ (k^ + k, k^) ~ (k^, kf,) for 
|k| < A. In all the results presented in this work we have taken A = 1 GeV 
and checked that this choice is in a region with a very weak dependence on A. 



4-2 Dependence of the Gluon Green's Function on the External Momenta 

After describing the technicalities of our method of solution we are now ready 
to show the behaviour of the NLL gluon Green's function. In Fig. 5 we have 
calculated f{ka, ki,, A) fixing k^ at 30 GeV and varying ka for the choice of 
parameter A = 3 and A = 5. We start by noting the complete agreement 
between the iterative and analytic results at LL. It is also interesting to notice 
how the angular integrated NLL gluon Green's function evolves from being 
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Fig. 4. Dependence of the NLL solution for the gluon Green's function on the 
parameter A for ka = 25 GeV, kf, = 30 GeV and A = 3. 

strongly peaked in the region ka — kh for small values of A to being more flat in 
ka for larger values. This behaviour is to be expected, since the BFKL equation 
can be reformulated as a differential equation in A, with 5*^^) (k^ — k^) as the 
boundary condition at A = 0. In agreement with this statement we obtain 
from Eq. (9) that / (k^, k;,, A = 0) = 5^^^ (k^ — k^). By comparing the plots in 
Fig. 5 wc also see how the gluon Green's function evolves from the boundary 
condition as A increases. 



Although the normahsation of the gluon Green's function is different between 
the LL and NLL cases, and indeed changes with A, the shape with respect 
to ka does not change too much in the ka — kh region. This behaviour is 
different when ka ^ k^, as can be seen in Fig. 6. In this region the gluon 
Green's function is much suppressed at NLL compared to the LL case. Similar 
results were obtained in Ref. [5] when the solution to the BFKL equation 
was studied using the Mellin transform of the angular averaged kernel, and 
a similar treatment of the running coupling terms as in our function ^ was 
considered, i.e. the f3o terms not rcsummcd into as- In Ref. [5] it was also 
shown that the resummation of these /3o terms improves the behaviour of the 
gluon Green's function. We will come back to this point in a forthcoming 
publication. 

In principle, it would be interesting to establish whether the Green's function, 
as obtained in the present approach, exhibits the exponentially suppressed, 
oscillatory behaviour predicted in Ref. [3] . For the accessible range of rapidities 
in our present numerical study it was not possible to verify this effect which 
would take place at large values of the ka/ ki, ratio and A. 

We now proceed to the study of the dependence on A in the next Section. 
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Fig. 5. ka dependence of the LL and NLL gluon Green's function at /x = fe^ = 30 GeV 
for two values of A. 

Dependence of the Gluon Green's Function on A 



In Fig. 7 we have plotted the evolution of /(ka, k^, A) with A for a specific 
choice of ka and /c^, both for the LL and NLL case. We see that for this choice 
of momenta, there is a significant difference in the A-dependence from about 
A = 3. With our choice of renormalisation scale of /i = /c^, the gluon Green's 
function rises slower with A at NLL than at LL. For the range in A studied 
in this paper no inestability in the growth of the gluon Green's function has 
been found. It is known from previous investigations in the literature [5] that a 
correct resummation of the 13q terms into the running of the coupling eliminates 
this possible problem. We will investigate this further in the future. 



We study a renormalisation scale dependence by, in Fig. 7, including the results 
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Fig. 6. ka dependence of the LL and NLL gluon Green's function at /i = fcf, = 30 GeV 
for A = 10 and a large range in the ka/kb ratio. 




Fig. 7. A-dependence of the NLL gluon Green's function evaluated for ka = 25 GeV 
and kb = 30 GeV. The central lines for both the LL and NLL result are obtained by 
choosing the renormalisation scale ^ = k^. The limits of the two upper (red) bands 
correspond to a choice of renormalisation scale /i = kb/2^ while for the bound of the 
two lower (yellow) bands, the renormalisation scale is /U = 2kb. 

for the choices /i = khjl and /x = 2kb. At LL, where the coupling does not run, 
we can still estimate a similar dependence to that at NLL by simply choosing 
the fixed value of the coupling at LL to be as(At = ^t), as(/i = kbl2) and 
as{n = 2kb). This generates the band around the LL result in Fig. 7. The 
upper (red) bands are obtained for /i < kf,, while the lower (yellow) bands 
correspond to fi > kf,- This renormalisation scale dependence is very big at 
LL and growing with A. For the particular selection of parameters, the scale 
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Fig. 8. The dependence of the gluon Green's function at NLL on the angle between 
k„ {ka = 25 GeV) and {h = 30 GeV) for the choice of A = 3 and A = 5. The 
renormahsation point is chosen at = kb. 



uncertainty is drastically reduced at NLL, with a slow growth with A. This last 
feature has its origin in that as A increases, more powers of are effectively 
resummed because there is more phase space available for emission, as was 
seen in Fig. 3. The reduction in the uncertainty at NLL with respect to the 
LL result shows how the predictive power of the theory has improved with 
the inclusion of radiative corrections. We will study this renormahsation scale 
dependence in a toy cross-section below, which will include an integration 
over a range of momenta. But before this, we first analyse in the next Section 
the angular behaviour of the gluon Green's function. 



4-4 Angular Dependence of the Gluon Green's Function 



One of the advantages of our method of solution is that we can perform studies 
of angular correlations in the NLL BFKL gluon Green's function. For example, 
in Fig. 8 we show the dependence of f{ka, kb, 9, A) defined in Eq. (13) for 
ka — 25 GeV and kj, — 30 GeV, and the two values A = 3 and A = 5. Again 
we see that for smaller A, the transverse momenta ka and kf, are strongly 
correlated, with the bulk of the contributions to the gluon Green's function 
coming from the region of small angles between the transverse momenta k^, k;,. 
We have included an investigation of the angular correlations of a toy cross 
section in Section 4.5. 
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Fig. 9. ^(A) calculated at LL and NLL, including the renormalisation scale depen- 
dence as detailed in the text. 

4.5 Toy Cross Section 

We now proceed to the study of the following quantity: 

^(A)= J ^ J ^/(k.,k,,A), (17) 

A:a>30GeV " fc5>30GeV * 

which, when the gluon Green's function is evaluated at LL, is proportional to 
the cross section at this accuracy. We therefore call it a "toy" cross section. 
A more complete study would require the use of the full NLL impact factors 
for different physical processes, a calculation which is out of the scope of this 
paper. Although the accuracy of the impact factors, LL, does not match that 
of the gluon Green's function, NLL, the behaviour of S{A) in Eq. (17) with 
A will still give an indication of the intercept we can expect at NLL. 

The evolution of the toy cross section with A is shown in Fig. 9 both for the 
LL and the NLL case. This plot is consistent with the fact that the NLL cor- 
rections decrease the intercept of cross sections. With the intention to estimate 
the renormalisation scale dependence in the LL curve, instead of keeping the 
strong coupling fixed at, for example its value at the lower integration limits 
as{n = 30 GeV) ^ 0.14, we fix it at Q;s(yU = kf,) as we do in the NLL cal- 
culation and in our study of the gluon Green's function in Sections 4.2- 4.4. 
Again we see, as in Fig. 7, that the renormalisation scale uncertainty is re- 
duced when going from the LL to the NLL curve for our particular choice of 
fi. When the renormalisation scale is chosen to be larger, the effective value of 
the coupling is reduced and the LL rise decreases, while the size of the NLL 
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Fig. 10. Average value of cos(0) for S{A) at LL and NLL as a function of A. At 
A = the two transverse vectors ka, kfe are completely correlated, whereas the 
BFKL evolution leads to a growing decorrelation with A, both at LL and NLL. 
This decorrelation is less at NLL than at LL. 

corrections, suppressed by one power of as, diminishes. The overall effect is 
that the difference between the LL and NLL curves is now smaller than for 
a lower choice of renormalisation scale. The initial decrease in S{A) in going 
beyond the LL approximation was already predicted in models adding running 
coupling effects to the LL evolution [8] . 

Finally, and in order to illustrate the feasibility of our calculations, in Fig. 10 
we have plotted the average value of cos(^^), with 9 being the angle between 
ka and k^, for the toy cross section as a function of A. We see that as A 
increases, and k^ become increasingly decorrelated, and we see that the 
effect is bigger at LL than at NLL. 



5 Conclusions 

In this study we have presented a first analysis of the behaviour of the NLL 
gluon Green's function as obtained from a numerical implementation of the 
method proposed in Ref. [6]. With the main purpose of showing the feasibility 
of our method we have taken a particular choice of renormalisation scale in 
the MS renormalisation scheme. Other choices and schemes can be considered, 
and work is in progress to study them. We have shown that the convergence 
properties are well understood, and in particular wc have demonstrated the 
infrared finiteness of our solution. We have also presented results on the evo- 
lution and angular dependence of the gluon Green's function and a toy cross 
section. The intercept obtained from this procedure decreases at NLL with 
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respect to the one obtained at LL, a trend which is in agreement with results 
in the hterature [4]. The magnitude of the change when going from LL to NLL 
depends on the choice of renormahsation scale and further studies are needed 
to draw stronger conclusions. 
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